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Abstract 

An approach which allows to include the corrections from non-orthogonality 
of electron states in contacts and quantum dots is developed. Comparison of 
the energy levels and charge distributions of electrons in ID quantum dot (QD) 
in equilibrium, obtained within orthogonal (OR) and non-orthogonal represen- 
tations (NOR), with the exact ones shows that the NOR provides a considerable 
improvement, for levels below the top of barrier. The approach is extended to 
non-equilibrium states. A derivation of the tunneling current through a single 
potential barrier is performed using equations of motion for correlation functions. 
A formula for transient current derived by means of the diagram technique for 
Hubbard operators is given for the problem of QD with strongly correlated elec- 
trons interacting with electrons in contacts. The non-orthogonality renormalizes 
the tunneling matrix elements and spectral weights of Green functions (GFs). 
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Introduction: In the description of tunneling processes tlirough quantum devices 
two approaches are useful: 1) the wave functions used for calculation of the Green 
functions fulfil the boundary conditions for the whole device, or 2) a subdivision of the 
system is made |]l], ^ ||, ^, ||, H and then the wave functions corresponding to different 
subsystems are, in general, not orthogonal to each other. The second approach is 
preferable if the strength of the interactions in different regions of the system differ 
considerably. In this framework tunneling arises due to the non-orthogonality of the 
wave functions from different subsystems. Prange [0 found in his investigation of 
SNS- and SIS-j unctions, that an overcomplete non-orthogonal basis set, allowing for 
the desirable separation leads to corrections from overlap integrals of the same order 
as the tunneling coefficients. Hence, it is essential to take the overlap between the 
orbitals into account. We also note the conclusion of Svidzinskii |^, that the tunneling 
Hamiltonian approach is useful only if one is interested in linear responses. 

We present results of a different approach based on the diagram technique for 
Hubbard operators within non-orthogonal basis f^. Generalized to non-equilibrium 
states, the method still allows for calculations in the language of model subsystems. 

Equilibrium: Consider a finite box with hard walls containing a barrier of finite 
height Vo see Fig. the Hamiltonian is 7i = + V . We approximate this system by 
two separate subsystems for which {0p,ep} and {(pq.Sq} are complete ON eigensystems 
of the 'left' (L) and 'right' (R) Hamiltonians Hl = ^ + Vl, p e L and Ur = + Vr, 
q E R corresponding to the potentials Vl = VQ{—x + b) + (Vq + V)Q{x — h) and 
Vr = iYo + V)Q{—x — a) + VQ{x + a) see Fig [1|; 9(x) is the Heaviside step function. 
Introduce the field operators ipii^) = J2pcr (^po-'Ppi^) i'Rix) = Hqa '^qo't'qi.^)] ^ 
denotes spin. The exact field operator is expanded as ip = ipA + i^B where ipA = 
i^L +'4'R ^iid is a reminder. Assuming ipB = yields the approximate Hamiltonian 

'Ha = T.ka^kalaO'ka + Epga [tpq^'la^^qa + H .c}j ,where tkk' = / '/'fc + ^) (l>k'^X] k = 
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p, q. Then tpq = OpqEq + Wpq, where Opq = J 0*0gdx defines the overlap matrix and 
Wpq = J 4>l(y — VR)(j)qdx, and similarly for the other matrix elements. Neglect the 
differences Wkk' whenever k, k' belong to the same contact. The operators Opo- are 
defined by ap^ = J2p' 0"^, J (f)*,{x)'ilj{x)dx + J2q' C>pq' 1 0*/(a;)-?/^(x)dx, and similarly for 
ttq^. Then, the anti-commutation relations are {a^a, ^1,^} = Of^,^,, where O^^, is element 
kk' of the inversed overlap matrix. Solutions of the Dyson equation g^^ = g^^ — 0~^W 
for O^^ = I (OR), I is the identity operator, O^^ 7^ / (NOR) compared to the 
exact solution are shown in Table |l[ As seen, the improvement achieved in NOR is 
considerable. However, in the proximity of the barrier height corrections from the 
reminder ipB should be taken into account when calculating the charge distribution, 
even though the energy levels estimated by NOR still are much better. 

Transient current: Transient current can be calculated from —e-^{NL), where 
A''^, is the number of carriers in the cylinder with top and bottom areas S and axis 
parallel to current. The tunneling matrix elements now contain the vector potential, 
so, in the Hamiltonian Ha substitute tkk' by Tkk' = (fc|(k — (e/c)A(t))^ + V\k'), the 
non-equilibrium tunneling coefficients. Whenever k, k' belong to the same contact we 
approximate Tky by its corresponding equilibrium value and neglect the differences 
Wkk', yielding for example Tpp/ Sp. 

Via the transformation api^q)^ = e~*'^-^(^)*^*)/po-(rgo-), which allows to introduce current 
states, we derive a tunneling current J ~ |Tp as a function of the applied voltage 
eV = fJ-L — fJ-R ^ ) ^ ^Jigj-g lY^Q phase ip = Lpi — Lp^. Putting O^y = Skk' when 

k, k' belong to the same contact, is consistent with the approximation. Note that with 
this assumption {Nl) = J2pa<Sp{npa), where Up^j = a^p^ap^ and Sp = Spp is a part 
of the overlap matrix, obtained by integration over the volume in which —e-j^{Ni) 
is calculated. The equations of motion for {up^) is -^{np^) = 2lmJ2qTpqe^'^^^^l^p^rqa) 
where (/J^r,,) = fpq{fp - /g) ^rrfe^' = ' ^p' = ^vq + ^m^i in the given 
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approximation and / is the Fermi function. The resuhing current formula in terms of 
densities of states Nl{^i) and Nr{^) then is J = V ■ A-ne^SS\f\'^NL{^i)NR{ii) = VR'^, 
(the factor 2 is due to spin). Thus, in this approximation the only change required is 
replacing T — T. 

Strong correlations in quantum dot: Physically, the coupling arises due to 
overlap of wave functions. When the system is close to the regime of Coulomb block- 
ade, the overlap is small and the interaction between subsystems is much weaker than 
the one inside the QD, i.e. the matrix element of tunneling |T| << = -Er^+i ^ E^^:, 
here i?r„ is an eigenvalue of the Hamiltonian of the QD, Hj\V) = Ey\T)] a = [r„, r„+i] 
is a fermion-like transition (a = [r„_|_i,r„]) which is described by Hubbard operator 
^r„r„+i ^ j5^a_ ^]-^jg gituation the states of the QD will be perturbed only slightly, 
therefore QD should be described by many-electron states. Here we will demonstrate 
how the technique developed in ref. [§] for correlated electron systems for thermody- 
namics can be extended to non-equilibrium phenomena. Since the Hamiltonian includes 
strong Coulomb interactions inside QD, the terms of the kind Vka,^i2fj.i^iiAa^'^iJ.2^tJ-i'^tJ-4 — 
Vka,fi2fisf^4{dli^2^f^'i^fj'i)'^^ka-^"' = T^a,a''^ka-^°'y ^^^o coutributc to the process of electron 
tunneling from the left contact to the QD and should not be decoupled in Hartree- 
Fock fashion. We include these interactions to the Hamiltonian of coupling (see ref. 
PI). The total Hamiltonian then is: 

^ = E ^kAaCka + E ErX^^ + E {TkaAaX"" + B.C.) 
kcr r ka,a 

where k = p,q [p & L and q E R). Any X-operator can be rewritten in terms of 
products of single-electron operators d, ct^ of QD and, therefore, using {cka, d)^} = Oj.^^^ 
one can show [H that {cka^X"-} = C'^J^,^((i^)''£:|"X^, where ^ is a Bose-like transition, 
and defines commutation relations between X-operators in the uncoupled system, 
{X^,X"-} = Strictly speaking, when the coupling is switched on, the many- 

electron states also become non-orthogonal to each other but the corrections contain 
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higher order products of C^cr,^- Hence, we neglect these corrections since we are only 
interested in first order with respect to transparancy. Following ref. we calculate the 
contribution to the current from 'left' electrons jj;!,} = —2eSSReYl,ka,a'^k<T*a^kcr ai^'^)^ 
where T^^^^.a — Tka,a + C^fca,/i('^M)*^a^-^r- The additional term comes from the anti- 
commutation of Cfccr and the QD Hamiltonian TId- Expressing the current in terms of 
retarded and 'lesser' GFs of the QD, G^^ and respectively, yields: 

4^ = 2e^5Im ^ {f!}:,0-,l^{d,YP'^ f^iek.) - \f!la?[GU^k.) + /L(e..)Gj,(u;)|,,j} 

ka,a 

For a transition a = [y,T] the expectation value = ({X^^,X^i'}) = + , 
i.e. it is a sum of the population numbers corresponding to the transition. They 
should be found from / Ga,a'\uj)duj. The system for Ga,a'\uj) is very cumbersome and 
therefore not given here. However, physics is seen from the form of the retarded GF: 
GU^) = .J-T-rZ)P^ ^ where r,(.;) = ri')(^) + rM(^), 7a(^) = li^u;) + 7^(0.), 
and the width r«(cj) = Y.T^^}*a9pA^)T^a,a- 9 is bare GF of 'left' electrons, t^'^cj) = 
T^^}*^gpa{uj)0~^^^{d^y'P^ and in T^^\ijj)., 'y^\uj) summation is over q,a. Thus, each 
single-electron intra-dot transition acquires width, which depends on the overlap of the 
wave function of conduction electron in the left (right) contact with energy near Fermi 
level fiiifJ'R) with those in-dot orbitals in transition a. 

In conclusion we have shown that, the separation of a device into auxiliary sub- 
systems unavoidably leads to eigenbases non-ortogonal to one another. This results 
in additional contributions to matrix elements of tunneling and in redistribution of 
spectral weights since part of charge is in 'intermediate' state. Precision of calculations 
is improved even in the most 'dangerous' region at the top of the barrier. 
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Exact 


OR 


NOR 


0.28 


0.012 


0.29 


1.06 


0.80 


1.08 


2.28 


2.07 


2.30 


4.04 


3.31 


4.18 



Table 1: Energy levels (mHartree) for barrier width 1 Bohr and height 4.1 mHartree. 
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Caption 

Fig. 1 The assumed potentials - 'left', 'exact' and 'right'. 
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Figure 1: Pransson et al. 
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